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ABSTRACT 

The Hipparcos data allow a major step forward in the research of 'moving groups' in 
the Solar neighbourhood, as the common motion of group members causes converging 
proper motions. Previous knowledge on these coherent structures in velocity space 
has always been limited by the availability, reliability, and accuracy of ground-based 
proper motion measurements. 

A refurbishment of Jones' convergent point method is presented which takes full 
advantage of the quality of the Hipparcos data. The original implementation of this 
method determines the maximum likelihood convergent point on a grid on the sky 
and simultaneously selects group members from a given set of stars with positions 
and proper motions. The refurbished procedure takes into account the full covari- 
ance matrix of the Hipparcos measurements instead of standard errors only allows 
for internal motions of the stars, and replaces the grid-based approach by a direct 
minimization. The method is tested on Monte Carlo simulations of moving groups, 
and applied to the Hyades. Despite the limited amount of data used by the convergent 
point method, the results for stars in and around the cluster-centre region agree very 
well with those of the recent comprehensive study by Perryman et al. 
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1 INTRODUCTION 

In the early years of this century, the acquisition of stel- 
lar positions and proper motions was one of the main aims 
of astronomical institutes. L. Boss' Preliminary General 
Catalogue (1910) revealed 'streams of stars' which showed 
converging proper motions (e.g., Kapteyn 1914, 1918; Ed- 
dington 1914). These streams were correctly interpreted 
as groups of stars with a common space motion: 'moving 
groups'. As projection effects are responsible for the converg- 
ing patterns, proper motions can be used to establish mem- 
bership of moving groups which are sufficiently extended 
on the sky. Therefore, convergent point studies have always 
been limited to nearby open clusters such as the Hyades 
(e.g., van Bueren 1952) and OB associations such as Scor- 
pius OB2 (e.g., Blaauw 1946; Bertiau 1958). Even though 
associations are gravitationally unbound, they form a coher- 
ent structure in velocity space: the initial velocity dispersion 
of association members is a few km s _1 (e.g., Mathieu 1986), 
similar to that of the gas in the parental molecular cloud. 
Due to expansion, tidal disruption by differential Galactic 
rotation, and stripping by passing molecular clouds, the size 
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of an association as well as its velocity dispersion grow with 
time, until the group eventually merges with the Galac- 
tic disk field star population. Associations generally have 
large physical dimensions without clear boundaries, and of- 
ten show signs of substructure in the form of subgroups (e.g., 
Blaauw 1964a), which may have different ages and kinemat- 
ics (e.g., Elmegreen & Lada 1977). 

In the pre-Hipparcos era, reliable high-accuracy proper 
motions were available for a limited set of stars only. In or- 
der to avoid zonal errors and other systematic effects, proper 
motion measurements were restricted to small samples of 
bright stars in fundamental catalogues or to small areas 
covered by a single photographic plate. Consequently, as- 
trometric membership of nearby associations, which cover 
tens to hundreds of square degrees on the sky, has always 
been limited to spectral types earlier than ~B5, thus hiding 
their low-mass stellar content. Consequently, many issues re- 
lated to the structure and dynamics of associations are still 
unknown (e.g., Blaauw 1964a, 1991). A well-defined mem- 
bership selection is the necessary first step to solve these 
open questions. Based on Hipparcos^] data, de Zeeuw et al. 
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The Hipparcos Catalogue (ESA 1997) is an all-sky astro- 
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(1999) have carried out a census of the stellar content of the 
nearby associations. Membership was established based on 
two independent selection methods. Both are based on the 
concept of common space motion: the 'Spaghetti method' 
uses proper motions and parallaxes to search for the com- 
mon motion of stars in velocity space (Hoogerwerf & Aguilar 
1999); the convergent point method is discussed in this pa- 
per. It is based on proper motion data only. 

Early algorithms implementing the convergen t p oint 

cf. 



2.2 



method are due to Charlier and Bohlin (1916; 
Smart 1968). Our version is basically a refurbishment of 
Jones' (1971) method (§2.2). Because of the limited com- 
puter power in the early 70's, Jones' original implementa- 
tion determined the maximum likelihood convergent point 
on a grid on the sky. We have replaced the grid-approach 
by a direct minimization routine. We also take into account 
the internal motions of the group members, as well as the 
full Hipparcos covariance matrix (§[2. 4). 

The convergent point method selects members based on 
proper motions and positions, but does not use any assump- 
tion about their distribution in configuration space. There- 
fore, the method can detect moving groups 'pur sang', the 
members of which fill a small volume in velocity space. How- 
ever, any attempt to find such moving groups in the full Hip- 
parcos Catalogue using proper motion data only is likely to 
yield a non-physical detection, as any convergent point has 
an associated 'moving group' covering large parts of the sky 
(§^^). In such cases, inclusion of parallax, radial velocity, 
and/or photometric data is essential to decide whether the 
space velocities of the candidate members are identical or 
not. Therefore, we have limited ourselves in this paper to 
test and apply the convergent point method to open clusters 
and associations. These two special types of moving groups 
(members of an open cluster or OB association have a com- 
mon space motion as well as configuration space position: 
they fill a restricted volume in phase space) allow the use of 
finite-size fields on the sky, resulting in an acceptable con- 
tamination by field stars (Q. Our simulated groups have 
30 pc diameters, requiring up to 20°x20° fields on the sky, 
and thus resemble moving groups/associations rather than 
classical open clusters. 

This paper is organized as follows. §^ describes the 
original and new implementation of the convergent point 
method. §^ deals with the optimization of the algorithm 
based on synthetic data, while ^ presents an application to 
the Hyades with Hipparcos data. §|E] discusses our results. 



2 CONVERGENT POINT METHOD 

2.1 Basic notation and definitions 

The Hipparcos Catalogue (ESA 1997) provides stellar po- 
sitions a, 8 and proper motions /i a cos<5, (is (in mas yr _1 ) 
in the equatorial coordinate system. Galactic coordinates 



metric catalogue containing absolute positions, proper motions, 
and parallaxes for 118,218 selected stars with accuracies of 
~1 mas (yr — 1 ). The limiting magnitude is V ~ 12.4 mag, and 
the completeness limit is V ~ 7.3—9.0 mag, depending on spec- 
tral type and position on the sky. The Catalogue also contains 
duplicity, photometric, and variability information. 



are denoted by t, b (ESA 1997, Vol. 1, §1.5). The trigono- 
metric parallax 7r (in mas) is related to distance D (in pc) 
according to: D = 1000/7T. The proper motions (i a cos8 
and fis are the components of the proper motion vector fj, 
which corresponds to the velocity tangential to the line of 
sight v tan (in km s _1 ) according to titan = A/i/ir, where 
A = 4.740470446 km yr s _1 is the ratio of one astronomical 
unit in kilometers and the number of seconds in one Julian 
year (ESA 1997, Vol. 1, table 1.2.2). For a given radial ve- 
locity « ra( j (in km s _1 ), the space velocity v (relative to the 
Sun, in km s _1 ) then follows from: 
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We denote by (X, Y, Z) the linear velocity components corre- 
sponding to v with respect to the usual rectangular equato- 
rial coordinate system in which the X-component is directed 
towards the vernal equinox, the V-component is directed to- 
wards the point on the equator with a = 90°, and the Z- 
component is directed towards the northern equatorial pole. 

In the absence of measurement errors and internal mo- 
tions (velocity dispersion), the proper motions of a set of 
stars with the same space motion v (a 'moving group') con- 
verge on the sky: 

, 7rusinA 

Ml = —j-, (2) 
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where v = \v\, and 

cos A = sin 8 sin S cp + cos 8 cos S cp cos(a cp 



(3) 



and A is the angular distance between the position (a, 8) of 
any star and the position (a cp , 8 cp ) of the convergent point. 
The convergent p oint denotes the direction of the projected 
space velocity v (§2.2). The radial velocities of the stars are 
given by u ra d = v cos A. 

By transforming the proper motion components [i a cos 8 
and (is into components /xy , directed parallel to the great cir- 
cle which joins any star and the convergent point, and /i±, 
directed perpendicular to the same great circle, the conver- 
gence of proper motions (eq. H) can be expressed as: 



with 

where 
tan# = 



nv sin X/A 
0, 
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sin 8 cos 8 
cos 8 sin 8 



jla cos 8 
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cos 6 tan S cp — sin 8 cos(a cp 



(4) 



(5) 



(6) 



Thus, the expectation value of /Uj. for any star is zero. 

The /ix-distribution contains information on the veloc- 
ity dispersion of the stars and on the measurement errors. 
The /iy-distribution contains, besides these two effects, a 
systematic perspective effect, depending on the angular dis- 
tance between the stars and the convergent point (sin A), 
and the individual parallaxes (n). Conversely, the assump- 
tion of a common space motion v for all stars allows the 
components jUii to be used to derive individual parallaxes 
(e.g., Jones 1971). A sophisticated method to obtain these 
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so-called secular parallaxes for a given set of moving group 
members is described by Dravins et al. (1997). 

Eq. (Q), and thus eq. (Q), is also valid if the stars are 
in a state of linear expansion (with expansion coefficient k 
in km s" 1 pc _1 ; e.g., Blaauw 1956, 1964b). In this case, 
however, the factor v in eq. (^) no longer solely corresponds 
to the space motion v of the stars, but also depends on the 
expansion velocity, i.e., v should be replaced by v' = v'(v, k). 
As a result, a kinematical detection of expansion requires 
radial velocities. It follows that if proper motions are used 
to determine the space motion or convergent point of the 
stars, one derives the genuine space motion plus the reflex 
of an expanding motion. 

2.2 Charlier's method 

Consider a set of stars with the same space motion v. For 
each star, the 3 equations of condition are (e.g., Smart 1968): 



- sin a cos a 

- cos a sin 8 — sin a sin S cos 8 
cos a cos 8 sin a cos 5 sin 8 




(7) 



The 3x3 matrix and the vector v (eq. |l|) follow directly from 
the observables a, S, n, [i a cos 8, /xa, v Ta ,d- The space motion 
components (X,Y,Z), the modulus v of the space motion 
vector v, and the coordinates (a cp ,<5 cp ) of the convergent 
point are related through: 



cos a cp cos 8 C p 
sin a cp cos <5 cp 
sin 8 cp 



(8) 



The parallax n can be eliminated from the first two lines of 
eq. (Q) by use of eq. (Q), which leads to the equation: 

aX + bY + cZ = 0, (9) 

where the coefficients a, b, and c depend on the observables 
a, 8, fi a cos 8, fis- The resulting equation of condition (g) can 
be solved for a set of stars in a least-squares sense to yield 
the ratios X : Y : Z, and thus the convergent point: 



arctan 



arctan 



VX 2 + Y 2 



(10) 



The third element of the space motion, its modulus v, can 
be determined only if radial velocities are available. This so- 
lution procedure is due to Charlier^] (1916). However, Char- 
lier's method incorrectly treats the least-squares coefficients 
a, b, and c, which depend on observables with measurement 
errors, as constants. This leads to systematic errors, which 
were sometimes 'corrected' for differentially (e.g., Seares 
1944, 1945; Petrie 1949a, b; Roman 1949). 

2.3 Jones' method 

Based on the work by Brown (1950), Jones (1971) devel- 
oped a maximum likelihood method for a simultaneous de- 
termination of convergent point and moving group mem- 

t Independently, Bohlin arrived at an equation of the form a' X + 
b'Y + c' Z = 0. However, Bohlin's and Charlier's equations are 
equivalent, except for different weights of the coefficients: a' = 
a/fi, b' = b/fi, and c' = c/fi. See Smart (1968). 



bership. This procedure avoids the conceptual problems of 
Charlier's method (§[2.2|). The basic idea is to determine the 
maximum likelihood coordinates (ct cp , S cp ) of the convergent 
point based on a comparison of the proper mo tion compo- 
nents with their expectation value of zero ({2.1). 



Jones starts with a sample of N stars with known posi- 
tions, proper motions, and corresponding errors. He overlays 
the sky a = 0°-360°, 8 = -90°-90° with a grid of trial conver- 
gent points (e.g., 24x12 cells of 15°xl5°), and numbers each 
grid point (e.g., i — 1, . . . , Ngrid = 288). Then, he follows 
the recipe outlined below. 

(i) Start at grid point i = 1. 

(ii) Assume that this grid point is the convergent point. 
The coordinates of this point are referred to as (a cp , 8 cp ). 

(iii) Calculate for each star j the error-weighted value t± 
of the component (eq. H): 



ti = 



in 



where to first order: 



(11) 



cr_L = (ae H || ) + (o> Q cos s cos 0) + (a M5 sin 8) , (12) 

and ae follows from eq. (0). Jones assumes that the quantity 
t± is distributed normally with zero mean and unit variance. 
Thus, the probability distribution for the given combination 
of fj,±j and <t_lj to occur is given by: 



Pj 



exp(-|tij)- 



(13) 



For a set of stars with the same space motion, each proper 
motion vector transforms into the component ; conse- 
quently, the expectation value for fj,± equals zero (eq. ^). 
(iv) Evaluate at the grid point the quantity X 2 : 



X A 



E 

3=1 



(14) 



As t± is distributed normally (cf. step [iii]), X 2 is distributed 
as x 2 with N — 2 degrees of freedom. 

(v) Determine X 2 for each grid point: repeat steps (i) 
through (iv) for i = 1, . . . , -/V gr id. 

(vi) The total probability P for the given set of calculated 
values of t± to occur is given by: 



N 
3=1 



(2tt 



.N/2 



exp(-iX 2 ). 



(15) 



Thus, the likelihood function C can be defined as: 

C= (27r) iv/2 P = exp(-iA 2 ), (16) 

and maximizing C is equivalent to minimizing X 2 . There- 
fore, the grid point with the lowest value of X 2 is defined as 
the convergent point (a cp ,<5 cp ). 

(vii) Given the number of degrees of freedom, evaluate 
the probability e that X 2 will exceed the observed value of 
X 2 by chance even for a correct model (e.g., Press et al. 
1992): 



r$[N-2]) Jxa 

where T(x) denotes the Gamma function for x > 0. 



(17) 
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(viii) If the computed value of X 2 is unacceptably high 
(e < e m in), reject the star with the highest value of \t±\, 
correct N^N — 1, and go to step (i). If the value of X 2 is 
acceptable (e > e m i n ), go to step (ix). 

(fx) The maximum likelihood convergent point is chosen 
as the grid point (a cp ,<5 cp ), and all non-r ejec ted stars in 
the sample are identified as members (cf. § |2.5| ). Thus, the 
determination of the moving group convergent point and 
membership are intricately linked. 

It is important to reject in advance from the sample 
all stars with insignificant proper motions (cf. Jones 1971). 
Insignificant in this respect means: 



t 



JL 



yj ill COS 2 5 + fj% 



< u 



(18) 



Because such stars have proper motions that carry 'no in- 
formation', the corresponding values for t± are small, inde- 
pendent of the position of the convergent point. Thus, these 
stars are likely never to be rejected. 

The convergent point method selects as members all 
stars with proper motion components fi± that are close to 
their expectation value of zero. However, not all stars with 
a small p,±_ necessarily are moving group members: stars 
at large distances generally have small proper motions, and 
correspondingly also small components fi±. Thus, the con- 
vergent po int membership selection is biased towards distant 
stars (cf. §ED). 



2.4 A refurbished convergent point method 

In view of the release of the Hipparcos data, we hav e modi- 

in three 



2.3 



fied the convergent point method described in 
ways, and extended it to include membership probabilities. 

2.4-1 Error propagation 

Unlike previous astrometric catalogues which quote (la) 
standard errors for the observables, the Hipparcos Catalogue 
provides the full covariance matrix for the measured astro- 
metric parameters. Thus, the propagation of errors (e.g., 
eq. |l2| ) must take the full covariance matrix into account. 
A general prescription to do so is given by ESA (1997, 
Vol. 1, §1.5). The specific case for the transformation of 
(fi a cos 5, (is) to (n\\,n±) is described in Appendix |aJ 

2.4-2 Internal motions 

With infinitely accurate measurements, moving group mem- 
bers would not necessarily have proper motions directed ex- 
actly towards the convergent point because of their velocity 
dispersion. Therefore, selecting stars with t± = will not 
identify all members. Consequently, when dealing with ac- 
curate proper motions, a certain amount of deviation of t± 
from should be allowed in the selection procedure. This is 
achieved by changing the definition of t± (eq. Illh to: 



t± = 



/'J 



v°±+< 



(19) 



where a* nt is an estimate of the one-dimensional velocity 
dispersion in the group, in proper-motion units. If (Ji nt is 



the one-dimensional velocity dispersion (in km s x ) and D 
is the distance of the group (in pc) one has: 



<T int = Aa* nt D/1000. 

The definition of t (eq. [ji]) is also modified: 
* - M 



(20) 



(21) 



Jones acknowledges that the assumption of a normal 
tx-distribution is critical to the use of the \ 2 distribution, 
but states that the rejection of individual stars depends 
on the much weaker assumption that pj decreases monotoni- 
cally as t± increases'. We have carried out Monte Carlo sim- 
ulations of moving groups with different distances, physical 
sizes, and velocity dispersions, in order to assess whether 
the observed ^-distributions are normal with zero mean 
and unit variance. Kolmogorov-Smirnov tests reveal no sig- 
nificant differences between the simulations and the model, 
independent of group distance, size, and velocity dispersion, 
though the inclusion of internal motions, i.e., usage of eq. 
( (lif ) instead of (fill), is essential. 



2.4-3 Direct minimization 



The original approach of evaluating X on a grid (§2.3) is 



limited in the sense that the accuracy of the position of the 
convergent point is restricted by the mesh size of the grid. 
Although a zoom-in strategy could overcome this problem, 
we can now drop the grid-based approach, i.e., improve the 
numerical implementation of Jones' method. In order to find 
the maximum likelihood convergent point, we apply a global 
direct minimization routine in two dimensions. This routine 
explicitly uses the analytically known derivatives dX 2 /da cp 
and dX 2 /dS cp of the objective function X 2 with respect 
to the free variables a cp and S cp (eq. |l4|). It returns the 
convergent point, as well as its uncertainties in the form of 
a 2x2 covariance matrix. 



2-4-4 Membership probabilities 

A natural membership probability should be defined in the 
/i||-/ij_-plane, as our membership selection takes place in this 
plane. For any star, the relevant observables are the proper 
motion components fin and p,±_ and the elements of the cor- 



responding 2x2 covariance matrix C (Appendix |Bj, eq. Bl), 
which describes the shape and orientation of the confidence 
region related to the Hipparcos measurements. A member- 
ship probability p can be defined as p = 1 — p CO nf, where 
Pconf is the minimum confidence limit of the confidence el- 
lipse that is consistent with p,±_ = 0. In Appendix [b], it is 
shown that p C onf = 1 — exp (— ^[fi±/a±] 2 ), so that the mem- 
bership probability p for a star is: 



exp( 



/'J 



(22) 



In view of the new definitions for t± and t (eqs [l9| and |2l| ) 
which take into account the velocity dispersion of the group, 
we modify the membership probability definition to: 



V 



ex P (-- 



2 



(23) 
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2.5 Discussion 

Whether a specific star is rejected or not depends, apart 
from (a cp , S cp ), not only on its individual membership prob- 
ability but also on the probabilities Pj of all other stars, 
since the product of the pj is considered (eqs |l^ and |l6|) . 
More specifically, if there are many stars with high pj (low 
tj_; very well converging proper motions) more low-pj (high- 
t±) stars can be tolerated. In other words: a moving group 
with a small velocity dispersion allows the method to include 
some non-members. However, this effect cannot be large as 
the velo city d ispersion of the group is accounted for in the 
model (§ |2.4.2[ ). 

Instead of using t± as membership indicator one could 
introduce a polar angle 4>t in the t-plane (tu = tsiacj>t, 
t± = t cos <j>t), and consider its distribution instead of f(t±) 
(see Appendix^] for details). Tests indicate that this alterna- 
tive membership determination procedu re g ives very similar 
results to the procedure presented in §§2.3-2.4. 



The current numerical implementation of the refur- 
bished convergent point method is not able to cope with 
samples containing two or more moving groups. The pro- 
cedure selects stars which are related to the lowest valley 
(measured in terms of X 2 ) in (a cp , 5 cp )-country. A detection 
of all minima, each one indicative of an over-abundance of 
high-pj stars, could simultaneously identify multiple moving 
groups. The statistical significance of each minimum could 
be assessed by considering the pj-distribution and compar- 
ing it with the expectation from field stars only. This method 
would be the two-dimensional analogue of Hoogerwerf & 
Aguilar's (1999) three-dimensional Spaghetti method. 



3 OPTIMIZATION OF THE ALGORITHM 

The convergent point method has several free parameters. In 
order to understa nd t heir importance, and determine their 
optimum values (§3.2), we apply the method to synthetic 
data of clusters sup erimposed on a kinematic model of the 
Galactic disk (§3.1). We have included the basic character- 
istics of the Hipparcos Catalogue in the synthetic data. 




300 240 180 120 60 
Galactic longitude / (deg) 

Figure 1. The effect of Galactic rotation and Solar motion on 
proper motion fi as function of Galactic longitude t, at 6 = 0°, 
for four different distances (D = 150,300,450, and 600 pc). The 
six vertical lines denote the specific choices for I of the cluster 
centres in the synthetic data. 



Table 1. Characteristics of the synthetic data. The radius of the 
spherical cluster is 15 pc, independent of its distance D. The total 
number of stars -/V mg + Nfi e \d in each field of view equals its size 
(Field size £xb) times 3. The latter number represents the mean 
stellar density per square degree in the Hipparcos Catalogue. 



D (pc) 


Field size ixb (deg 2 ) 
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Afield 


150 


20x20 


200 


1000 


300 


15x15 


100 


575 


450 


10x10 


50 


250 
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10x10 


25 


275 



3.1 Construction of synthetic data 

The synthetic data set consists of 24 distinct configurations: 
we vary the Galactic longitude I of the cluster centre from 
t = 0° to 300°, in steps of 60°; we assume a constant Galactic 
latitude b = 0°. We vary the distance of the cluster centre 
from D — 150 to 600 pc, in steps of 150 pc. We represent 
each synthetic cluster by JV mg stars distributed in a sphere 
of radius 15 pc with a constant volume density. As the clus- 
ter distance increases, the fictitious number of members ob- 
served by Hipparcos decreases due to the completeness limit 
of the Catalogue. Because we consider clusters with a fixed 
intrinsic size, we decrease the size of the field of view with 
increasing cluster distance. The number of field stars Afi e id 
is chosen such that the total number of stars in the field is 
consistent with the Hipparcos Catalogue mean stellar den- 
sity of 3 stars per square degree. The size of the field of view 
as well as the number of cluster members iV mE have been 
chosen based on results obtained for the nearby associations 



by de Zeeuw et al. (1999) . Table |l] summarizes the properties 
of the different configurations. 

3.1.1 Cluster kinematics 

The measured proper motion of each cluster star can be de- 
composed into (1) a systemic streaming motion of the cluster 
with respect to its own standard of rest, (2) a component 
due to velocity dispersion, (3) the reflex of the Solar motion 
with respect to the Local Standard of Rest (LSR), and (4) 
the effect of differential Galactic rotation. First, we take ef- 
fects (l)-(3) and observational errors into account; then we 
add Galactic rotation to the proper motion. 

Ad (1): the streaming velocity of the cluster is randomly 
drawn from a constant-density sphere in velocity space with 
a radius of 10 km s _1 . Ad (2): the velocity dispersion is 
assumed to be isotropic, and each of the three independent 
components is drawn from a Gaussian distribution with a = 
2.0 km s _1 . Ad (3): we adopt for the Solar motion with 
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150 - 



s 100 
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log(e) 




log(e) 



log(e) 



Figure 2. Results of the membership selection procedure applied to synthetic data with a cluster centred at (£,b) = (300°, 0°) at a 
distance D = 150, 300, 450, and 600 pc. The panels refer to the different cluster distances as indicated. Each panel shows the number of 
accepted cluster stars (solid) and field stars (dashed) as function of loge (eq. [n]). The dashed vertical line at the right-hand side of each 
panel denotes e m i n = 0.954 or loge m i n = —0.020. Note the different vertical scales of the panels. 



respect to the LSR the Hipparcos value of Dehnen & Binney 
(1998): (U, V, W)q = (10.00, 5.23, 7.17) km s -1 . 

We translate the synthetic three-dimensional velocities 
and positions into observables, after which observational er- 
rors are added: the positions are assumed to remain un- 
changed, while we disturb the parallax and proper motion 
by amounts which are randomly drawn from Gaussian distri- 
butions with a — 1 mas (yr _1 ), consistent with the median 
errors of the Hipparcos data. Then, we add Galactic rotation 
to the proper motion (item 4) using the first-order formu- 
lae (e.g., Smart 1968). We adopt for the Oort constants the 
Hipparcos values of Feast & Whitelock (1997): A = 14.82, 
B = -12.37 km s" 1 kpc" 1 . Finally, the covariance matrix 
elements for the synthetic member are directly copied from 
a star which is randomly drawn from the subset of stars in 
the Hipparcos Catalogue which have positive parallaxes. 



in terms of the Schwarzschild velocity ellipsoid (e.g., Binney 
& Merrifield 1998). This model is simple but adequate for 
our purposes. Based on Hipparcos data, Dehnen & Binney 
(1998) measured the first two moments of the velocity dis- 
tribution in the Solar neighbourhood as function of B — V 
c olour, as well as the Solar motion with respect to the LSR 
(§3.1.1). We use their results to characterize the shape and 



orientation of the velocity ellipsoid. 

We randomly draw the required number of field stars 
from the subset of stars in the Hipparcos Catalogue which 
have coordinates in the corresponding field of view as well 
as positive parallaxes. We retain all Hipparcos measure- 
ments, but replace the proper motions by synthetic values 
which are projections^] of space velocities randomly drawn 
from Schwarzschild's model, given the corresponding B — V 



3.1.2 Field kinematics 

The synthetic field star distribution is consistent with the 
characteristics of the Hipparcos Catalogue as well as with 
the kinematical properties of the local Galactic disk. We de- 
scribe the velocity distribution in the Solar neighbourhood 



§ The projection requires 7r > 0. Although other, more compli- 
cated, methods could be used which do not require tt > 0, it is 
not a priori clear that such efforts give more realistic results, es- 
pecially because only 4 per cent of all Hipparcos stars has 7r < 
(or no parallax at all). 



© 0000 RAS, MNRAS 000, 000-000 



Moving groups in the Hipparcos Catalogue 7 



0.8 



0. 




0.3 



o. ; 



0.6 



0.8 




Figure 3. The fractions / mg (solid) and / ne id (dashed) of ac- 
cepted cluster and field stars, respectively, as function of e for 
the specific configuration with cluster centre coordinates (£, b) = 
(300°, 0°) and distance D = 150 pc. The value e min = 0.954 is 
indicated by the dashed vertical line. The two lowest / mg -lines 
correspond to the two 'observations' outside the central band in 
Figure ||. 



colour. Galac tic ro tation is added to the proper motions as 



described in 83.1.1 



Following this procedure, the characteristics of the syn- 
thetic field stars are consistent with the characteristics of the 
Hipparcos Catalogue (e.g., the V magnitude, parallax, and 
spectral type distributions) as well as with the kinematics 
of the Galactic disk in the corresponding field of view. De 
Zeeuw et al. (1999; their §§3.4.3-3.4.5) extensively discuss 
Dehnen & Binney's results and the procedure described here 
in relation to similar Monte Carlo simulations. 



3.2 Application to synthetic data 

We have constructed 20 random realizations of a synthetic 
cluster superimposed on a field star distribution for each of 

and Table |l|) . For each of 



3.1 



the 24 specific configurations (§ 
these 480 synthetic data sets, the cluster convergent point 
and members were determined using the refurbished con- 
vergent point method (§2.4). After extensive experiments, 
we choose to fix t mln = 1.7 (eqs [t| and |l|), and to use 
for ai n t (eq. ^) the value which was actually used in the 
construction of the synthetic data (2.0 km s" 1 ). The lat- 
ter choice is motivated by the fact that the results of the 
membership selection turned o ut t o be rather insensitive to 
the precise value of a lnt (cf. §4T). Our choice for i m i n is 
motivated by the median distance of the associations dis- 
cussed by de Zeeuw et al. (1999) of ~300 pc. This value 
combined with (Ji nt = 2.0 km s _1 gives a* t = 1.4 mas yr _1 . 
For c M = 1 mas yr _1 , eqs (|l8|) and ( plf ) with t mln = 1.7 
imply that all stars with /i < 3 mas yr~ <ire rejected, which 
is identical to rejecting all stars for which the proper motion 
is less than the familiar 3a. The remaining free parameter 



Figure 4. Filled squares indicate the differences between the po- 
sitions of the real and determined convergent points for the 20 
Monte Carlo realizations with (£, b) = (300°, 0°) and D = 150 
pc. The dashed ellipses denote the 3cr u ncertainty regions of 
the determined convergent points (§ 2.4,3| ), The contours indi- 
cate constant values of X 2 , linearly spaced by AX 2 = 2000.00 
with a minimum value X 2 = 2.19 at (Ao cp , A<5 cp ) = (0°, 0°) or 
(acpAp) = (73?897, -1?747), for the 'perfect cluster' with Solar 
motion and Galactic rotation but without streaming motion. The 
two 'observations' outside the central band correspond to the two 
lowest / m g-lines in Figure jj| 



£mm (step [viii] in §2.3) represents the stop criterion in the 
selection process. We determine its optimum value based on 
Monte Carlo simulations. 



3.2.1 Effect of I and D 

The convergent point method kinematically detects clusters 
using proper motions only. A detection requires the group to 
stand out from the field star distribution sufficiently clearly. 
Field stars are background and foreground stars as well as 
non-cluster members at approximately the cluster distance. 

The distinction between background and foreground 
stars and cluster members is simple for nearby clusters. 
This is illustrated in Figure [l], which shows the relation be- 
tween Galactic longitude £ (at b — 0°) and proper motion 
due to the distance-independent effect of Galactic rotation 
(Feast & Whitelock 1997) and the distance-dependent reflex 
of the Solar motion (Dehnen & Binney 1998) for an object 
at D = 150, 300, 450, and 600 pc. The figure shows that (1) 
for given £ the measured proper motion of a star is an indi- 
cator of its distance, and (2) this indicator is more sensitive 
at smaller distances, especially near £ — 150° and 300°. 

The cluster streaming motion enables to distinguish be- 
tween field stars at the cluster distance and genuine clus- 
ter members . In the simulations, this streaming motion 
(item 1 in §3.1.1) has a maximum value of 10 km s - . 
isotropic case, the expectation value for the tan- 



In the 

gential streaming velocity equals 157r/ 



5.89 km s" 
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(jf — 3.14 . . .), which translates to an amplitude in proper 
motion of fi — 5.89 it /A = 1.24 tx mas yr -1 ; this gives 
fi ~ 8.3 mas yr™ 1 for D = 150 pc, but only /x ~ 2.1 mas yr -1 
for D — 600 pc. Besides the magnitude of the streaming 
motion, its orientation is also important. A 'mis-aligned' 
streaming motion could cause the proper motions of the 
group members to vanish, or could cause them to mix in- 
distinguishably with the Galactic disk distribution, reducing 
the astrometric cluster detection probability. 



3.2.2 Effect of e 

All results of the membership analysis and convergent point 
determination for the 24 specific configurations can be un- 
derstood completely in terms of the arguments presented 
above. As an example, Figure shows the results for the 
20 synthetic data sets with the cluster centred at {I, b) = 
(300°, 0°) for D = 150, 300, 450, and 600 pc. The figure shows 
the selected number of field stars (20 dashed lines) and clus- 
ter stars (20 solid lines) as function of loge, the stop param- 
eter. As e increases, the cluster detection significance level 
increases (eq. ^), which means that more stars have been 
rejected in order to reach the lower X 2 . Ideally, this rejec- 
tion concerns field stars only. This would lead to solid lines 
which are constant at the number of cluster members down 
to e = 1, and to decreasing dashed lines, vanishing at e — 1. 

Figure |^ shows that one generally starts with fewer clus- 
ter members than iV mg because some of them are rejected 
due to their insignificant proper motions (eqs |l^ and pl[ ). 
Furthermore, the number of selected field stars does not 
vanish at e = 1. This indicates that the membership lists 
are contaminated by field stars ('interlopers'). These stars 
have proper motions which are consistent with the projected 
space motion of the cluster. Elimination of interlopers re- 
quires radial velocities, photometric data, and/or distance 
information. The number of selected cluster stars is not con- 
stant with e but decreases due to the erroneous rejection of 
members, caused by, e.g., velocity dispersion, observational 
errors, and/or an unfavorable streaming velocity. Especially 
for large cluster distances, many members have sometimes 
been rejected before an acceptable solution is reached: in 
these cases, the combined reflex of the streaming motion, 
Galactic rotation, and Solar motion is insufficient to kine- 
matically detect the cluster using proper motions only. 

Figure | shows, for (£, b) = (300°, 0°) and D = 150 pc, 
the fractions / mg and /field of selected cluster and field stars, 
respectively, as function of e. The convergent point method 
confirms ~80 per cent of the cluster members, but with a 
contamination of ~20 per cent of the total number of field 
stars in the field of view. Simulations show that these per- 
centages are practically independent of the total number of 
cluster or field stars in the simulation. 

The optimum value for the stop parameter e, e m in, is the 
value which returns the largest number of cluster members 
with the smallest contamination by field stars. It follows 
from Figures ^ and |^ that e m in should be chosen close to 1, 
independent of cluster distance. This conclusion is supported 
by results of the other synthetic data sets. We have chosen 
e min = 0.954 (loge min = -0.020). 



3.2.3 Convergent point 

Figure ^ compares the positions of the real convergent point 
of the cluster and the one determined by the convergent 
point method for the 20 Monte Carlo realizations with 
(£,b) = (300°, 0°) and D = 150 pc; we define Aa cp = 
c*c P ,obs — a?cp and AS cp = 5 cp ,obs — 8 cp . The 'observations' 
form an elongated structure, extending over several tens of 
degrees, in the Aa cp -A5 C p-plane. 

This effect is explained by the contours in Figure ^, 
which represent constant values of X 2 (eq. |l^). The con- 
tours correspond to the 'perfect cluster' at (£, b) = (300°, 0°) 
and D = 150 pc: 200 stars with proper motions due to 
the reflex of Solar motion and Galactic rotation only. Each 
of these proper motion vectors defines a great circle. Be- 
cause the proper motions for all stars are nearly parallel, the 
corresponding great circles are nearly parallel. As a result, 
the proper motions define quite accurately a great circle on 
which the convergent point must be located; however, the 
precise location of the convergent point on this great circle 
is rather uncertain. The contours in Figure ^ show exactly 
this: the convergent point has a large freedom, in terms of 
X 2 , to move along a great circle whereas excursions in the 
direction perpendicular to this great circle are not allowed 
(e.g., Bertiau 1958; Perryman et al. 1998). 



4 APPLICATION 

4.1 Hyades open cluster 

Perryman et al. (1998; P98) combined Hipparcos proper mo- 
tions and parallaxes with ground-based radial velocities in 
order to establish membership for the Hyades open cluster. 
Their sample contains all stars in the Hipparcos Catalogue 
in the field 2 h 15 m < a < 6 h 5 m and -2° < S < 35° with 
■k > 10 mas. P98 define as members all stars for which the 
space velocity coincides with the Hyades centre-of-mass mo- 
tion to within the 99.73 per cent confidence region ('3<r'). 
For stars lacking radial velocity information, the decision is 
based on tangential velocity only. P98 start with 142 secure 
pre-Hipparcos members in the central region of the Hyades. 
Convergence of their procedure was achieved after two it- 
erations, resulting in 218 members, 179 of which were pre- 
viously known (all with known radial velocities), and 39 of 
which are new (18 of these have a known radial velocity). 
The P98 member selection is generous; only few genuine 
Hyades members have probably not been selected, whereas 
some field stars are likely to be present in the member list. 

Figure 8(a) in P98 shows the positional distribution of 
the 218 members. P98 consider four distinct components (r 
is defined as the three-dimensional distance to the cluster 
centre): (1) a spherical 'core' with a core radius r c ~ 2.7 pc 
and half-mass radius r\ ~ 5.7 pc; (2) a 'corona' extending 
out to the tidal radius r t ~ 10 pc; (3) a 'halo' consisting of 
stars with r t Si r ^ 2r t which are still dynamically bound to 
the cluster; and (4) a 'moving group population' of stars with 
r <; 2r t which have similar kinematics as the central parts. 
The halo and the moving group population are explained by 
evaporation and diffusion of core and corona members due 
to the effects of the Galactic tidal field and stellar encounters 



(§4.2). The core and corona contain 134 stars (r < 10 pc), 
the halo contains 45 stars (10 < r < 20 pc), and the moving 
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Table 2. The number of selected Hyades members by P98 (Perryman et al. 1998), B99 (this study; cr int = 2.0 km s^ 1 ), B99 and P98 
(B99 I P98; I for intersection), P98 and not B99 (P98 M B99; M for minus), and B99 and not P98 (B99 M P98), for four different 
regions around the Hyades cluster centre (0 < r [All], r < 10 pc ['core' + 'corona'], 10 < r < 20 pc ['halo'], and r > 20 pc ['moving 
group population']). The comparison between B99 and P98 is made for three different sets of P98 members: '3cr' members are all stars 
for which the space velocity and the Hyades centre-of-mass motion agree to within the 99.73 per cent confidence region; for '2<r' members 
these quantities agree to within 95.4 per cent, while for 'lcr' members they agree to within 68.3 per cent. 
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Figure 5. The number of Hyades members found by B99 (this study; solid black), by P98 (dashed black), the intersection B99 I P98 
(dotted black), the difference P98 M B99 (dashed grey), and the difference B99 M P98 (solid grey), as function of the (one-dimensional) 
velocity dispersion <7i nt used in the selection procedure. From the top row down: 'lcr' P98 results (confidence limit p con f = 0.683), '2cr' 
P98 results (p C onf = 0.954), and '3tr' P98 results (p C onf = 0.9973). From the left column right: stars with r < 10 pc ('core' + 'corona'), 
stars with 10 < r < 20 pc ('halo'), and stars with r > 20 pc ('moving group population'). 
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Figure 6. Results of the membership selection of the Hyades by P98 ('3<r' members; top row), B99 (this study; CT; nt = 2.0 km s 1 ; 
second row), B99 and P98 (B99 I P98; third row), P98 and not B99 (P98 M B99; fourth row), and B99 and not P98 (B99 M P98; fifth 
row). From the left column right: (1) positions on the sky; (2) vector point diagram (fi a cos8,fig); (3) distance histogram (D = 1000/7T; 
D in pc, 7r in mas) with number of stars indicated; and (4) membership statistic or probability as function of three-dimensional distance 
r (pc) from the Hyades cluster centre: p (0 < p < 1; eq. versus r for B99 and B99 M P98; the dimensionless quantity c, defined by 
c = z T S _1 z, where T denotes taking the transpose, z is the difference vector between the observed space velocity of the star and the 
Hyades centre-of-mass motion, and S is the sum of the two corresponding covariance matrices, for P98 and P98 M B99. In case z is 
i/-dimensional, c is distributed as x 2 with v degrees of freedom: for v = 2 (no radial velocity information), c = 11.83 corresponds to a 
confidence limit p CO nf — 

0.9973 ('3cr'), and for v = 3, c = 14.16 relates to p conf = 0.9973 ('3ct'). 



group population contains 39 stars (r > 20 pc). Kinematic 
assignment of membership beyond ~2rt is difficult, and thus 
less reliable. 

Figure 8(b) in P98 shows the velocity distribution of 
the 218 members. The expected one-dimensional velocity 



dispersion is 0.2-0.4 km s~ , too low to be resolved: the 
observed velocities are fully consistent with parallel space 
motions of the members, with a (one-dimensional) velocity 
dispersion of ~0.3 km s _1 due to a combination of internal 
motions, undetected binaries, and measurement errors. 
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We apply the convergent point method to the Hyades 
with D — 45 pc, t m in = 1.7, and e m i n = 0.954. Figure ^ 
shows the results as function of the (one-dimensional) ve- 
locity dispersion cri nt used in the selection. It follows from 
Figure [B| that: (1) in order to reproduce the P98 results, 
a non-zero value for o"j nt is required; for <7i nt J; 1.5 km s _1 , 
the numbers of selected stars are comparable to the numbers 
found by P98; (2) despite the fact that P98 have made use of 
additional parallax and radial velocity data, the convergent 
point method is able to retrieve nearly all P98 members in 
the inner 20 pc of the cluster centre, without introducing 
many spurious members at the same time. Most of the stars 
that are not picked up by the convergent point method have 
a low P98 membership probability (cf. Figure ^|) . The con- 
vergent point method picks up a number of field stars in the 
outer regions of the cluster, r <; 2r t ~ 20 pc (cf. Figure [(J). 
Parallaxes and/or radial velocities are required in order to 
reject these stars as candidate Hyades members. 

Table ^| gives the comparison of this study (B99) and 
P98 results for cr int = 2.0 km s" 1 . Fi gure q shows the results 
of the membership selection and the comparison with the 
'3<r' P98 results for <7i n t = 2.0 km s _1 . The convergent point 
method correctly identifies 203 of the 218 P98 Hyades mem- 
bers. Only 15 P98 Hyades members are not picked up by the 
convergent point method: 4 of the 15 lie between 9.30 and 
16.47 pc from the cluster centre, but have a low membership 
probability. The same is true for the remaining 11 stars, at 
r > 26.70 pc. Furthermore, only 4 of them have a known ra- 
dial velocity. Eighty-seven new candidate Hyades members 
are identified by the convergent point method. Five of them 
have r < 10 pc, but rather low membership probabilities 
(p < 0.60). Seven stars have 10 < r < 20 pc. The remaining 
75 candidates have r > 20 pc, combined with low mem- 
bership probabilities as well as predominantly small proper 
motions and parallaxes (D ^ 50 pc), and are distributed 
uniformly over the sky. This indicates that these stars are 
most likely field stars, which can be rejected as candidate 
members using parallax and/or radial velocity information. 

We find X 2 = 245.574 (0.853 per degree of freedom), 
and (a cp ,S cp ) = (95?54 ± 0?43, 7?28 ± 0? 19) with a cor- 
relatio n coe fficient p — —0.79 on the measurement er- 
rors (§ [2. 4, 3 ). Projecting the Hyades centre-of-mass motion 



(table 3 in P98) gives corresponding convergent points of 
(a cp ,S cp ) = (97?91,6?66) for the 134 stars with r < 10 pc, 
and (97? 96, 6? 61) for the 179 stars^J] with r < 20 pc. When 
we apply the convergent point method to the 218 P98 mem- 
bers as starting set, we end up with 213 members having 
X 2 = 151.712 (0.702 per degree of freedom), (a cp ,8 cp ) = 
(97?81 ± 0?52, 6?74 ± 0?21), and p = -0.84. The 5 rejected 
stars have r > 30 pc, and 4 of them have no radial velocity 
information. The different convergent points seem to show 
a large spread. However, it can be explained fully by the 
errors: the confidence ellipses are elongated along the great 
circle which connects the convergent points and the Hyades 
cluster centre, and the spread in the convergent points is ex- 
actly along this great circle. Accordingly, the position of the 



^ We neglect the small difference in cluster centre coordinates 
between the 134 stars within 10 pc and the 180 stars within 20 pc 
of the cluster centre (table 3 in P98), and correspondingly only 
have 179 stars within 20 pc of the cluster centre. 
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Figure 7. The 3392 'comoving' Hyades members with n > 
10 mas: sky distribution (bottom), parallax distribution (mid- 
dle), and membership probability distribution (top). The two 
solid dots in the bottom panel denote the convergent point 
(acpAp) = (95? 54, 7?28) and its mirror point (275? 54, -7? 28). 



convergent point is well-defined in the direction perpendic- 
ular to this g reat c ircle, whereas it can easily shift along the 
great circle (§ 3.2.3). This is illustrated by the rejection of the 



last star in the application of the convergent point method to 
the 218 P98 members; this step shifts the convergent point 
from (a cp ,<5cp) = (96?70,6?91) to (97? 81, 6? 74). 



4.2 Hyades supercluster 

Although it has often been claimed that the Hyades open 
cluster is related to a kinematically similar 'supercluster' 
(e.g., Hertzsprung 1909; Eggen 1992; Chereul, Creze & Bi- 
enayme 1998), these results are not commonly accepted (e.g, 
Ratnatunga 1988; Skuljan, Cottrell & Hearnshaw 1997). 
This extended moving group population is generally ex- 
plained by the evaporation of the Hyades open cluster 
through (1) relaxation by internal stellar encounters, (2) the 
Galactic tidal field, and (3) passing interstellar clouds (e.g., 
Terlevich 1987). Although all these effects are indisputably 
taking place, it is not to be expected that all escaped mem- 
bers have velocities consistent with the cluster space motion: 
stars beyond the cluster tidal radius are subject to Galactic 
tidal forces and to the shearing effect of differential Galac- 
tic rotation, and therefore can have systematically distorted 
motions. A selection criterion different from velocity, e.g., 
angular momentum, might therefore be more appropriate to 
identify escaped members (e.g., Innanen, Harris & Webbink 
1983; Helmi, Zhao & de Zeeuw 1998; P98). In short, estab- 
lishing the significance of the reality and the interpretation 
of superclusters and all-sky stellar streams is an important 
issue requiring a thorough analysis, beyond the scope of this 
paper. 

Nonetheless, an interesting question is: 'is there a signif- 
icant over-abundance of low-tx stars for the Hyades conver- 
gent point over the whole sky?' We have applied the conver- 
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gent point method to the entire Hipparcos Catalogue, sub- 
ject to the constraint n > 10 mas, to search for 'comoving' 
Hyades members. Out of 22951 targets, we find 3392 such 
stars, the overwhelming majority of which (94 per cent) has 
r > 20 pc. However, their distribution on the sky is com- 
pletely random, their parallax distribution peaks sharply 
towards ir — 10 mas, and their membership probabilities 
are distributed uniformly (Figure Q). All these facts are in- 
dicative of a non-physical grouping. Indeed, searching for 
all-sky comoving stars at random cluster velocities (conver- 
gent points) yields identical results. An all-sky search using 
the ACT/TRC catalogues (§^) is practically feasible, but 
is prone to give even less well-interpretable results because 
these catalogues (1) each contain 10 times more stars, (2) 
have proper motion accuracies which are 2-3 times worse 
than those of Hipparcos, and (3) do not allow the use of 
parallaxes to make the tt > 10 mas cutoff (giving in practice 
a 50 times larger sample). Obviously, proper motion data 
alone is insufficient for all-sky moving group analyses; an 
in-depth investigation of multi-band photometric, spectro- 
scopic, and/or parallax data is required to bring down the 
number of interlopers significantly. 



5 DISCUSSION 

The Hipparcos Catalogue (ESA 1997) has opened a new 
era of stellar kinematical research. The high quality of the 
data, its homogeneity, and the absence of systematic errors 
down to the 0.1 mas level, make the Hipparcos proper mo- 
tions ideally suited for a reassessment of the existence and 
membership of moving groups in the Solar neighbourhood. 
As Hipparcos parallaxes do not yield useful constraints on 
individual tangential space velocities for stars at distances 
D <; 350-500 pc, a common motion of distant stars is best 
detected by converging proper motions. The nature of the 
Hipparcos data, specifically the full covariance matrices pro- 
vided in the Catalogue and its high quality which sometimes 
allows to marginally resolve internal motions, combined with 
increased computer power over the last decades, have neces- 
sitated a refurbishment of the convergent point method that 
was originally developed by Jones in 1971. This method is a 
maximum likelihood procedure which detects moving groups 
and selects the corresponding members from a given set of 
stars with positions and proper motions. The method selects 
all stars which have proper motions that are consistent with 
the maximum likelihood convergent point, and assigns indi- 
vidual membership probabilities. The method is known to 
work well for nearby groups, such as the Hyades, which span 
a large enough area on the sky so as to allow the proper mo- 
tions to show their converging pattern (see Cudworth 1998 
for a review) . The method is biased towards distant stars, as 
these have generally small proper motions. Therefore, such 
stars are more likely not to be rejected as members. 

Extensive Monte Carlo simulations have been used to 
fix free parameters. These simulations show that typically 
~80 per cent of all cluster members can be retrieved, com- 
bined with a contamination of field stars which amounts to 
~20 per cent of the total number of (field) stars in the sample 
(Figure |^). We presented an application of our method to 
the Hyades with Hipparcos positions and proper motions. 
The results for stars near the centre of the cluster agree 



with those of Perryman et al. (1998), who used parallaxes 
and radial velocities as well to determine membership. An 
a posteriori inspection of the Hipparcos parallaxes of stars 
that are selected by the convergent point method but not 
by Perryman et al. confirms that nearly all of them lie at 
large distances from the Hyades cluster centre (Figure [|). 
The majority of these stars have low to medium member- 
ship probabilities, consistent with them being field stars. 

The abovementioned bias can also be removed by com- 
bining the convergent point method with other selection pro- 
cedures, such as Hoogerwerf & Aguilar's (1999) 'Spaghetti 
method' which also uses parallax information. De Zeeuw et 
al. (1999) have done so, and successfully applied the com- 
bined procedure to nearby OB associations with Hipparcos 
data. As expected, the combination of the two methods effi- 
ciently eliminates the biases in the individual methods (fig- 
ure 4 in de Zeeuw et al.). De Zeeuw et al. used a two-step 
strategy to determine membership, thereby increasing the 
sensitivity and power of the convergent point method. In 
the first step, they identified the association using the early- 
type stars in the Catalogue, with the projected space mo- 
tion, i.e., convergent point, as result. As the number of field 
stars among the early-type stars is generally several factors 
smaller th an the num bers assumed in our Monte Carlo simu- 
lations (§§ ^.2.l| - ^2~3[ Table 0), the convergent point method 
was able to detect associations out to distances of ~650 pc. 
In the second step, additional association members of later 
spectral type were selected from the remaining stars in the 
Catalogue by assuming the convergent point of the late-type 
members is equal to that of the early-type members. 

Several parametric and non-parametric methods for 
finding moving groups have recently been developed (Chen 
et al. 1997; Figueras et al. 1997; Chereul et al. 1998; Hooger- 
werf & Aguilar 1999; see also Dehnen 1998). All of them use 
Hipparcos positions, proper motions, and parallaxes; some 
of them also require ground-based radial velocities. The lat- 
ter methods are currently of limited use, as a uniform set of 
homogeneous radial velocities for the majority of stars in the 
Hipparcos Catalogue is not available (cf. Udry et al. 1997); 
the lack of reliable radial velocities for early-type stars is 
especially severe, so that these methods cannot be used for 
detection and membership selection of OB associations. 

The positions for ~10 stars with V ^ 11 mag in the 
Astrographic Catalogue have recently been combined with 
those in the Tycho Catalogue. The resulting ACT and TRC 
catalogues list proper motions on the Hipparcos reference 
frame with accuracies ~3 mas yr _1 (Urban et al. 1998; H0g 
et al. 1998). The ACT/TRC catalogues are presently not 
the ideal catalogues to apply the convergent point method 
to: in order to bring down the number of interlopers to an 
acceptable level, a large amount of additional spectros copi c 
and/or multi-band photometric data is indispensable (§4.2). 
Future astrometric satellite missions such as GAIA will pro- 
vide micro-arcsecond astrometry and accurate multi-epoch 
multi-band photometry down to V ~ 20 mag as well as ra- 
dial velocities accurate to ~3 km s~ down to V ~ 15 mag 
for ~10 9 stars (e.g., Perryman, Lindegren & Turon 1997; 
Gilmore et al. 1998). Such data will allow the detection of 
common motion of stars throughout the Galaxy! 
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APPENDIX A: FROM (/x a cos 8, /j, s ) TO (/x h /x x ) 

For a given convergent point (a cp ,S cp ), the transformation 
of (/i a cos <5 = ii a »,fL&) to can be carried out using 

the general recipe outlined in ESA (1997, Vol. 1, §1.5.2). 
The transformation uses a 4-dimensional vector a: 



a 
6 

Pa* 
(IS 



( «' 

8 



(Al) 



where a* = a cos 8. The covariance matrix of a, Cs., is given 
by Ca = J/ C a J/, where T denotes taking the transpose, 
J / denotes the Jacobian matrix of the transformation f, and 
C a is the covariance matrix of a. Using eqs (|B|) and (^|), we 
find that J/ = d(a* , 8, fin, p±)/d(a* , 8, fi a *, is given by: 
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V 
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cos f 
sin ( 



° \ 




(A2) 
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sin (a C p — Oi, 

ffsinO I 1 

\da* J 

(tan <5 tan 5 C p -f 



(tan 8 — tan S cp cos[a cp — a]) , 



cos 8 sin e 
sin(a C p — a) 
cos[a cp — a]) . 



(A3) 



APPENDIX B: MEMBERSHIP PROBABILITY 

The Hipparcos Catalogue provides the five astrometric pa- 
rameters a = (a cos 8, 8, n, fi a cos 8, /i«) T , where T denotes 
taking the transpose, together with the 5x5 covariance ma- 
trix C. Appendix [a] describes how the proper motion com- 
ponents iia cos 8 and lis can be transformed into /ny and fj,±, 
with the corresponding 2x2 covariance matrix: 



a I 
peril cm 



„2 



(Bl) 



where p denotes the correlation coefficient between and 
H±. We denote the 'measured' proper motion components p,^ 
and by TIjf and JI±. The confidence ellipse with confidence 
limit p C onf in the py-px-plane is defined by the locus of 
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points (/Ltii, which solve the equation (Press et al. 1992, 
p. 690-693): 

~W\ H±--PI)C- 1 ( ^1 ~5L )=A, (B2) 
v 7 y fix - Mi / 

where the parameter A determines the size of the ellipse, 
and is related to p con f according to: 



Pconf 



A/2 



dt t 



-1+I//2 



exp(-t), 



rW2) J:! 

with ^ the number of dimensions. After definition of: 

— mii — T^TT _„j rr-_Mx — mx 



and 



Mi 



(B3) 



(B4) 



I 1/2 



eq. (B2) can be written as: 

m]| 2 — 2pfb\\jTi_ + ji± 2 = A(l - 

from which follows that: 

Mx = MM|| ± [(l - P 2 ) (A - / 

A natural membership probability p can be related to the 
confidence limit p con f corresponding to the smallest confi- 
dence ellipse that is consistent with mx = in the point 
(M||,Mx) — (m||>Mj_)- Solving the equations: 



(B5) 
(B6) 



Mx' = -Mx/fx, 

L IU M||=M||' 



(B7) 




a exp (§ ffl2 ) 



(C2) 



where erf (a:) is the error function (e.g., Abramowitz & Ste- 
gun 1974, eq. 7.1.1) and a = ly sin^t. For t\\ > (< 0), the 
function f(<f>t) peaks at <f>t = 5 (%")j the width of the peak 
increases with decreasing i.e., with increasing cluster 
distance. Thus, an alternative to minimizing X 2 (eq. |l^) is 
minimizing: 



X,i, 



N 
3=1 



!J t,3 



N 
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3=1 



(C3) 
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where m(i an d Mj_ are related through eq. (B5), gives the 
solution (/ij,,Mx) — (±pVA~, ±\/A), with the result that: 

2 



A = 



(B8) 



The confidence limit p CO nf is given by eq. (B3) with v = 2: 
Pconf = 1 — exp (— | A). Thus, for a given star with mea- 
sured proper motion components Mil and Mx and covariance 
matrix C (eq. Bl) the confidence limit p con f of the smallest 



confidence ellipse that is consistent with fi± = is given by 
Pconf = 1 — exp( — i [/ix/crx] 2 ). The membership probability 
p for this star is then given by (cf. eq. E3): 



Pconf = exp ( 



(B9) 



APPENDIX C: ALTERNATIVE PROCEDURE 

Instead of using t± as membership indicator under the as- 
sumption that its distribution is normal with zero mean and 
unit variance, one can also work in t -space, where <f>t is the 
polar angle in the f-plane: tu = t sin <j> t and t± = tcos(j> t . If 
t± and t|| are distributed normally (§ |2.4.2 ) with unit vari- 



ance and mean zero and i\\ , respectively, it follows that: 



d^tdt 

+{tsin t - i||} 2 ], (CI) 



f^t,t)d<t> t dt = -X^texp(-i[{tcos^ t } 2 + 

Z7T 



such that: 

= f dt/(^,t)d^ = -^exp(-II|| 2 ) 
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